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CL|| McKean and Trubowitz [2] showed that the theory of the KdV equation 

d d'^ dg 

^ : at^^^'^^ = -Q^9{x,t)-%g{x,t) — {x,t). 
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is intimately related to the geometry of a related hyperelliptic curve of infinite 
\0 • genus, the Bloch spectrum Bg_^ of the operator 

o ■ 



o 



where gt = g{x,t). As was known classically, Bg^ is independent of t, when g{x,t) 
evolves according to the KdV equation. Our purpose in this paper is to develop a 
g3 ! theory of finite difference operators and their Bloch spectra and isospectral flows 
^1 which mimics the KdV theory. The basic idea of this paper is to use the theory of 
6jO' the periodic Toda chain of length A^. Here again, the periodic Toda chain can be 
• understood in terms of a finite genus hyperelliptic curve and isospectral deforma- 
>■ . tions, as van Moerbeke discovered. For instance, see [3]. So one would like to see 
^ I what the relation of the Toda hierarchy is to the KdV hierarchy, how the conserved 
^ [ quantities are related and so forth. A start on these matters has been obtained 
by Toda in [4]. In this paper, the idea is that if we choose the initial data for 
the periodic Toda chain very carefully, then the evolution of this data under the 
various equations of the Toda hierarchy looks similar to the evolution of / under 
the KdV hierarchy. Given /, we will find a canonical choice of the initial data of the 
Toda equations so that the flow of this initial data under the Toda hierarchy looks 
like the flow of / under the KdV hierarchy, at least to high accuracy. This choice 
will be given by an asymptotic series in A^~^. The main result of this paper is the 
formulation and partial verification of the Conjecture given below. The method 
studied here also gives a way of producing analytically approximate solutions to 
the Toda chain hierarchy, at least conjecturally. My motivation is to use this case 
as a model for constructing finite genus models for the KP hierarchy, which I have 
studied in [1]. I also hope to use these methods to develop discrete models of the 
sine-Gordon equation, the non-linear Schrodinger equation and other infinite di- 
mensional integrable systems that are related to isospectral flows. At the end of 
this paper, I discuss using the methods developed in this paper to construct a finite 
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difference equation to numerically solve the KdV equation. This finite difference 
equation has many conserved quantities. 

We first review the KdV hierarchy from an isospectral point of view. Let T be 
the translation operator defined on real valued functions on R by 

T{g){x) = g{x + l). 

Let C be the set of functions g so that g is analytic on R and T{g) = g, i.e. g is 
periodic with period 1. If g & C, then we can define an operator 

LgW = -i^" + g{x)il;{x) 

on the space of infinitely differentiable functions i/j on R. We define the Bloch 
spectrum Bg of Lg to be the set of (A, a) G C x C* so that there is a non-zero 
function i/j with Lg{'i(j) = Xt(j and T{'i(j) = aif). The KdV hierarchy is a sequence of 
non-linear differential operators -Di : C C so that for any i if we have an analytic 
function g{x,t) and we define gtix) to be g{x, t) and 

, , d 

Di[gt) = -^gt, 

then Bg^. is independent of t. 

The Toda hierarchy is quite similar to the KdV hierarchy. Let us fix an positive 
integer N. If A : Z ^ C, then we can again define the translation operator T by 
T{A){n) = A{n + N). Let Cn be the set of (A, B) with T{A) = A and T{B) = B. 
Given S) e Civ, we can form an operator on the space of all functions ip : Z ^ C 

by 

L^A,B)W{n) = {-iP{n + 1) + A{n)iP{n) + B{n)iP{n - 1))N\ 

We define the Bloch spectrum B(^a,b) of to be the set of (A, a) e C x C* 

so that there is a non-zero function with L(^A,B)i'^) — A'i/' and T{tl)) = ai/j. Then 
there is a sequence of non-linear difference operators Di,k{A, B) and D2.k{A, B) so 
that if At and Bt are in Cjv and At{n) and Bt{n) are differentiable functions of t 
and we have 

and 

D2,k{AuBt) = ^^Bu 

then B(^At,Bt) is independent of t. 

The Toda hierarchy can be defined in the following way: We first inductively 
define a sequence of complex valued functions di on Z as follows: We take di = 
for z < and do = 1. We also take di{0) = and finally 

di{n + 1) = di{n) + A{n)di-i{n) + B{n)di-2{n). 

The first Toda equation is defined by 



THE TODA HIERARCHY AND THE KDV HIERARCHY 



3 



and 

D2,i{A, B){n) = B{n){A{n) - A{n - 1))N. 

To define the k*'^ member of the Toda hierarchy {k > 1), we define a sequence of 
functions a^-i, . . . , ao by descending induction on p by 

fc-i 

ap{n) = dk-p{n + k) - dk-p{n) - ^ ar{n)dr-p{n + r). 

r=p+l 

We define 

Di,k{A,B){n) = {ai{n)B{n+l) - B{n)ai{n - 1))N 

and 

D2,k{AB){n) = B{n){ao{n) - ao{n - 1))N. 

The numbers di{N) are conserved quantities of the flow. If necessary, we denote 
di{N)hjdi{N,A,B).' 

Given g E C, our question is whether it is possible to find N and {A, B) e Cn so 
that the Toda flows with initial data A and B look like the KdV flows for g and so 
that Bg looks like Ba,b- Rather than deflne precisely what we mean by looks like, 
we will just try to solve the problem through a series of attempts. We first define 

1 

We will use this notation consistently throughout the paper. We first consider the 
problem of making Ba,b look like Bg. We try the following: 

A{n) = 2 + e^g{en) 

■"^ While it is not necessary for the rest of the paper, here is the algebro-geometric interpretation 
of the Toda lattice equations above. The curve Ba.b is a double cover of the A line. It can be 
compactified by adding two points P and Q over A = oo and A^^ = 2; is a parameter at these 
two points. Further, the function a has a pole of order N dX P and a zero of order N aX Q. 
Let h{z) = az^ , regarded as a function of z in a neighborhood of P. To simplify notation, let's 
assume that /i(0) = 1. If A and B are generic with /i(0) = 1, then given a generic effective divisor 
D of degree g, the genus of Ba,Bj we can find unique functions Sn in L(D + n(P — Q)) so that 
the top coefficient of the Laurent expansion of s„ is one. Thus sq = 1 and s{N) = a. We define 
di{n) to be the z*'^ Taylor coefficient of z^Sn at P. Thus di{N) is the i*'^ Taylor coefHcient of h{z) 
expanded as a power series in z. In particular, di(N) does not depend on D, but only on the curve 
Ba,b- It is an easy consequence of Riemann-Roch that Sn, s„-|-i, Sn—i, and Xsn satisfy a linear 
relation, which turns out to be of the form 

As„ + Sn+l = AoSn + BoSn-l 

where A^) and Bjy are functions which depend only on the linear equivalence class of D. Further, 
there is a divisor Dq so that A^^^ is the original A and Bu,, is the original B. It is easy to find 
a recursive formula for the di{n) in terms of the Au and Bd- We can map the curve Ba,b to 
the Jacobian Jg of Ba,b by j{x) = Dq + P — x. The derivative of j at P is an element of the 
tangent space of Jg at Dq, which can be extended to a translation invariant vector field on Jg. 
The first Toda equation is just the derivatives of the functions Ajj and Bd as functions on Jg with 
respect to this translation invariant vector field. The higher Toda equations come in a similar way. 
Instead of using the tangent vector to j{BA,B), one uses vectors in the higher osculating subspaces 
of j{BA,B) at P. The calculation of the Toda equations from this point of view is essentially worked 
out in [1,§3]. The Toda equations are the case M = 1 of that paper. The paper shows how to 
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and 

B{n) = -1. 

Suppose that Lg^il^) — A'0 and Ti^ip) = a^ip so that (A, a) G Bg. Define 4>{n) = ipine). 
Plugging this A and B into 

-(/)(n + 1) + A(n)(/)(n) + B{n)ct){n - 1) 
^(A,B)(0)(n) = 

and using Taylor's theorem to expand '4>{e{n + k)) around en, we see that 

(f 

LA,B{(l>){n) = --^'^i^'^) + 9{^n)ilj{en) + 0(e). 
Since Lg{^) — Xijj, we have 

LA,B{<f>)in) = ^<f>{n) + 0{e). 



and T{(f)) = acf). Thus (f) is an approximate eigenfunction of La,b and T. It is at least 
plausible that there is a true eigenfunction of La,b and T with LA,B{<Pe) — Xe(pe 
and T{(f)e) = cte^e and 

lim (Ae, CKg) = (A, a). 

Conversely, if A is fixed and N is large, it is reasonable to assume that if we have 
(A, a) G i3iv, then there are points of Bg nearby. Thus thus Bg looks like Ba,b with 
this choice of {A,B). More precisely, we conjecture that if K is any constant and 

SK = {{X,a):\X\<K}, 

we can make the two sets Sk H Bg and Sk H Ba,b be as close together as we want 
by making large. 

The problem with this approach is that the above choice of A and B is somewhat 
arbitrary. If we let 

A{n) = 2 + e^{^g{en) + q{en)) 

and 

B{n) = -1 + e\-g{en) - q{en)), 

the same argument shows that Bg looks like Ba,b with this choice of {A, B). How- 
ever, the first conserved quantity 

N-l 1 

d^{N) = V A{k) ^2N + e / -g{x) + q{x) dx + 0{e^) 
fc=o 2 

depends on our choice of q. That is, we know what A + Bisto order e^, but we do 
not know how to define A and B. But the conserved quantities depend on A and 
B, not just on A + B. So given g, we want to have a canonical choice of A and B 
in terms of g. Then we can hope that the conserved quantities of the A and B will 
just depend on g. However, we want this canonical choice of A and B to be in some 

; ; J- ] — ; rp„ i„ n 
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We will first concentrate on the first Toda equation: 

L»i,i(A5)(n) = {B{n) - B{n+1))N 

and 

D2,i{A, B){n) = B{n){A{n) - A{n - 1)N. 
Let us consider two functions ft{x) — f{x,t) and ht{x) — h{x,t). Set 

At{n) = 2 + e'^{f{en,t) + h{en,t)) 

and 

Bt{n) = -1 + e^i fien, t) - h{en, t)) 

and assume that f{x, 0) = g{x)/2. With this choice of A and B, Eg still looks like 
the Bloch spectrum attached to Aq and Bq. In general, i32/t looks like BAt,Bf Let's 
consider the equations 

^^A,in) = D,,,iAB)in) + Oie') 

and 

^^Bt{n) = D2,i{A,B){n) + Oie'). 

That is we ask that At{n) and Bt{n) be approximate solutions to the Toda flow 
problem. This means that BAt,Bt presumably do not move much as t changes 
and hence 62/^ also does not change much with t. Plugging in and using Taylor's 
theorem, we get 

j^{f + h){en) = -^{f-h){en) + 0{e) 

and 

±(f-h)ien) = -^{f + h)ien) + Oie). 

So if / and h are going to produce solutions to the Toda flow that are even approx- 
imately correct, then we should ask that 

l(f + h){x,t) = ^^{f-h){x,t)+0{e) 

and 

±(f-h)ix,t) = -^{f + h){x,t) + 0{e). 

Obviously, we cannot achieve this by taking f = h, which corresponds to the case of 
B = —1. However, we can take h = and then we do get an approximate solution 
to the Toda flow with the choice 

(1) At{n)^2 + e\f{en,t)) 

and 
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with 

Equations (1) and (2) define a canonical choice for A and B in terms of / = g/2 
modulo e'^. Of course, this is a rather uninteresting solution. 

(1) and (2) only define A and B to order e^. We next ask if we can modify (1) 
and (2) to be canonical to order e^. 

(4) At{n) = 2 + e^fien, t) - e^^3{f){en, t) 
and 

(5) Bt{n) = -1 + e2/(en, t) + e^^sif){en, t), 

where $3(/) is an unknown polynomial in / and the derivatives of /. We want to 
choose $3(/) in such a way that 

j^At^{Bin)-Bin + l))N + 0{e^) 

and 

j^Bt = B{n){A{n) - A{n - 1))N + 0{e^) 

for some choices at least of f{x,t). It is not obvious that such a $3 should exist. 
Using Taylor's theorem to evaluate A{n + i) and B{n + i), we should have 

(Q) ^ _ ^Mlle = + e(-^ - ^^^) 

^ ^ dt dt ^ dx 2dx^ dx ' 

and 

df d^sif) ^_4f d'f d^sif) . 
^ ' dt dt ^ dx 4da;2 dx 

In particular, adding (6) and (7), we see that f{x,t) should satisfy 
(8) jJ{x,t) = -^f{x,t) + 0{^). 



Now since 



we have 



|/(x,*) = -A;(,,t) + 0(.), 



<i<l>8(/)_ <<*3(/) 



dt dx 
by using the chain rule to evaluate 

dt 

So assuming that (8) holds, then (6) will hold if 

d$3(/) 1 d-^f 
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Thus we can take 
Thus if (8) holds and 



At{n) = 2 + e\f{en,t))-e\-\^{en,t)) 



and 



B,[n)^-l + e\f{en,t)) + e\-\^{en,t)), 



then A and B satisfy the Toda equations to order e^. Next we attempt to continue 
this process. That is, we try to find a function 

$4(/) 

of / and the derivatives of / so that if we define 

(9) At{n) = 2 + e\f{en,t)) - e\--J-{en,t)) - e'Mf)M 
and 

(10) S,(n) = -l + e2(/(en,t)) + e3(-^£(en,t)) + e4$4(/)(en,t), 

then A and B satisfy the Toda equations to order e^, provided that f{x,t) satisfies 
some condition. To find this condition on /, we plug (9) and (10) into 

(11) j^At = (Bin) - B{n + 1))N + 0{e') 
and 

(12) j^Bt = B{n){A{n) - A{n - 1))N + O(e^). 
Adding (11) and (12), we see that / must satisfy: 

Assuming (13) and using the fact that 



dt dx 

some computation shows that 



+ C(e), 



dMf) _ Kdf 
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So we see that if 

*4(/) = -^/^ 

then both (11) and (12) wiU hold if (13) holds. 

When we say that / satisfies (13), we mean that (13) is satisfied in the asymptotic 
sense. More generally, consider the equation 

(14) ^ = + evE'i(/) + . . . + e-^M) + 0{e^+'), 
at ax 

where the are polynomials in / and the x-derivatives of /. Of course, it can be 
extremely difficult to solve 

(15) |.,|+,W,(/) + ... + e"*„(/). 

So instead we try an asymptotic solution of the form 

(16) fix, t) = foix, t) + h{x, t)e + ... + Ux, t)e^. 

To solve our problem, substitute (15) into (16) and multiply everything out and 
denote the coefficient of e'' by Vk- Then we can easily find fk for k from to n so 
that Vk = for k from to n. Indeed, all we have to do is to solve 

where the are polynomials in the known fo,... , fk-i and their derivatives. 
Given initial conditions, we can then easily solve these equations. So an asymptotic 
solution in this sense to (13) will solve (11) and (12) to order for t in any given 
bounded region. That is, for any given T, there is an integer N and a constant C 
so that if n > AT, then 

\^^At-{B{n)-B{n + l))N\<Ce^ 

and 

\^Bt - B{n){A{n) - A{n - 1))N\ < Ce^ 

if \t\ < T. Of coiusc, these estimates may break down over long time intervals, i.e. 
C may depend on T. However, this kind of asymptotic analysis provides a tool for 
guessing the behavior of the equations of the Toda hierarchy, although rigorously 
connecting the asymptotic analysis and the behavior of the equations is probably 
extremely difficult. 

I conjecture that this process of finding can be continued indefinitely. At any 
rate suppose we define 

, / 1 d^fjx) , f(^\2^ii^^4 

7680 dx^ ^64^^''^ dx >^ 
_ J_ d^fjx) 3 d-'fixf 1 df{x)d^f{x) 5 
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and / satisfies 

df _ df{x) ^ ^ 1 d^f{x) ^ 1 ^^^^df{x)^^2 
dt dx 24 dx^ 2 dx 

.( 1 d^fjx) . 1 .2df{x) 1 df{x)dy{x) 1 d^fjx) 4 

1920 ^16^^^ rfx 32 dx dx^ '^48^^^ dx^ ' 

J_lf{xf 1 df{x)d'f{x) , 
^ 128 dx^ 128 da;3 ^ 

. 1 _ 1 d^fix)d^fix) 

^ 322560 dx'^ 128 da; 1536 da;^ dx^ 

^ 384-^^ ^ dx^ ^64 dx ^ dx2 ^64-^^ ^ dx 
1 dV(a:) 1 dfix)d^f{x) 6 

3840'^^^^ da;5 1536 dx dx"^ 

If we define 

(17) At{n) = 2 + e2(/(en, t)) - e^R{f){en, t) 
and 

(18) St(n) = -1 + e\f{en, t)) + e^R{f){en, t), 
then 

= -(S(n + 1) - B{n))N + 0{e^) 

j^Bt = B{n){A{n) - A{n - 1))N + 0(e9). 

These computations were carried out using Maple, as is described in the Appendix. 

Next let's consider the second pair of equations in the Toda hierarchy. Actually, 
the answer is more comprehensible if we try to solve the equations: 

(19) (i^i,2 + 2D^,^){At, Bt) = j^At + 0(6^) 
and 

(20) (^2,2 + 2^2,1) (Ai, Si) = As^ + o(e9). 

Now define A^ and B^ by the equations (17) and (18). In order for (19) and (20) 
to hold, we must have that 

e'^ = (Di,2 + 2L»i,i + L»2,2 + 2D2,i)(^t, Bt) + 0{e^). 
Plugging in our definition of A and B and using Taylor's theorem many times, we 
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dt ^ 4 dx^ ' dx ' 

3 .<lm__lff ,2df{x) 1 d^fjx) 
dx^ S^^""^ dx 64 dx^ 
^ ll df{x)dyix) ^^^ 
16 dx dx'^ 
(21) ^ ^ 3 d'^fjxf ^ 3 dfjx) d^f{x) ^^^ 

64 64 dx"^ 

1 d'fjx) _ ]_dl{xf 49 d^fjx) d^fjx) 

2560 rfx^ 64 dx 768 rfx2 rfx^ 
5 .2d^f{x) 7 4f(a;) .d''f{x) 9 .grf/la;) 



64"'^ ^ c?a;3 32 cZx ^ ^ dx'^ 32"'^ ^ da; 

il ff^.fJM + 35 d/(a;)dV(a;) e 
640-^^^ dx^ 768 cZa;4 ^ 

Conversely, if (21) is valid and we use (17) and (18) to define A and B, then (19) 
and (20) are valid. 

Next we repeat this process for the third set of equations in the Toda hierarchy. 
Here the result is that if 

df .5 .d^f{x) 15 .2df{x) 1 d^f{x) 



dt 4 dx^ 2 dx 16 dx^ 

^ 5 df{x) dV(g;) 4 
2 (ix dx'^ 

1 d'fjx) _ b^xf 155 d^fjx) d'fjx) 
192 dx'^ 4 dx 192 dx2 dx^ 
_ 45 .2^!M _ ^dm d^ 15 3#(^) 
32-^^^^ dx3 16 dx ' dx^ ^ 8-'^^ dx 
U dV(x) 47 d/(x) dV(x) e 
64 ■^^^^ 96 dx dx^ ''^ 

and we use (17) and (18) to define A and B, then 

(19) (i^i,3 - 2i?i,i + 2i?i,2)(At, Bt) = j^At + 0(e9) 



and 

(20.) (^2,3 - 2r>2,i + 2D2,2)(^t, Bt) = j^Bt + 0{e^) 

Finally, for the fourth Toda equation, we have that 

(19) (Di,4 + 4L>i,i - 2L»i,2 + 2Di,3)(^t, Bt) = j^At + 0{e^) 

and 

d 
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if 

^ = ((_LflM. _ ^^dfjxf 35 d^fjx) d^fjx) 
dt 64 dx^ 8 dx 16 dx^ dx^ 

35 .2^!/M 35 (//(x) 35 34f(x) 



8 (ix"^ 2 dx dx"^ 2 

^ fi^im , 21 df{x)d^ fix) , 

16''^ ' dx^ 16 dx dx^ ' 
This suggests the following conjecture: 

Conjecture. There are polynomials o,i^d 'ifkj{f) in f and the derivatives 

of f so that if 

L 

Rif) = J2^k{f)e'' 



and 

and we define 
and 



L+l 

^.•(/) = E*m(/)^' 

Atin) = 2 + e\f{en, t)) - e^R{f){en, t) 



Bt{n) = 2 + e\f{en, t)) + e^R{f){en, t) 
and f satisfies the equation 

| = Z,(/) + 0(e^+i), 

then At{n) and Bt{n) satisfy the j*^ Toda equations to order e^^^. Further, by 
taking suitable linear combinations of the Zj{f), we can produce asymptotic series 
whose leading terms in e are the KdV hierarchy if L is large enough. 

For the rest of the paper, let us assume this conjecture is true. Then we can also 
calculate out the conserved quantities asymptotically. First, let's consider di{N), 
the first conserved quantity. We have 

JV-l 

d,{N)=J2Mn). 

n=0 

By the Euler-Maclaurin summation formula, we have that asymptotically 

d^{N) = N{ ['{2 + e^fix) - e'Rif)) dx). 
Jo 

So we can write out the first few terms of the expansion of di (N) : 

di(N)^-(2 + e f f(x)dx+-e^ f f(xfdx + —e^f f(xfdx + .... 
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In particular, the following quantity is conserved asymptotically: 

Cl= [ f{x)dx + \t^ I f{xfdx+^t^ I f{xfdx + .... 
Jo ° Jo 32 Jo 

Next consider d2{N). We can write 

AT-l N-1 N-1 



n=0 n=0 n=0 



Again the Euler-Maclaurin summation formula allows us to calculate each of the 
terms as an integral. In order to make the formulas look nicer, we consider to 
conserved quantity 



(22) C2 = -l{d,{N) - I + f ) + ^-j^Cl) - \CI^ 



Then we have an expansion for C2 



f{xf dx + e'{-J^' f{xf dx-^l^' fix) dx) + o(6^). 



Thus if we solve a given Toda equation asymptotically to high enough order, we 
can expect 



f{xf dx + e\\ f{xf dx-^j^' fix) 



to be conserved to order e^. 

The third conserved quantity is C3 which has the expression: 

dsiN) - —(-16 + 72e - SGe^ - 24Cle + GOCle^ 

-36Cle3 + 12Cl2e2 + 18C2e^ + 2Cl^e^ + 18C2€^ - 9ClC2e^). 
We have the following expression for C3 

Let's use (19) and (20) to try to find a numerical method for solving KdV. We 
first define a(n) = N'^iAin) - 2) and 6(n) = iV2(S(n) + 1). Then we define 

La,bik) = + 1) - 26(A;) - a(A; + 1) + a(A; - 1) 

and 

Ma,bik) = 2a(/c) - 2a(/c - 1) - 6(/c + 1) + 6(/c - 1). 
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and 

Ga,b{k) = -2b{k)a{k)+2b{k)a{k - 1) +a{k)^ -a{k-l)^ + bik)b{k + l) - b{k)b{k - 1) 

+ {-b{k)a{kf + b{k)a{k - lf)/N^. 

(19) and (20) become 

(23) ^ik) = N{LaAk) + e'Fa,b{k)) 
and 

(24) ^{k) = N{Ma,b{k) + e'Ga,b{k)). 

Let's suppose we have a function f{x) we want to use as initial data for KdV where 
/ is periodic with period 1. Define c(n) = f{n/N). Then we can mimic (17) and 
(18) to define a(n) and b{n) to any desired degree of accuracy in e. Thus let's define 

A(6)W = ivW'^ + D-K'^-i)), 

So we can approximate the first derivative of / by 

1 



A(c) - ^e2A3(c). 



For instance, we could set 



(25) a{n) = c{n) + -/{/^{c){n) - \^l^\c)) + \e'c{nf - ^^'^\c). 
and 

(26) b{n) = c{n) - -/{l^{c){n) - \e''^'{c)) - jcinf + ^e'^^c)) 

and use these a and b as initial conditions for (23) and (24) to get functions a{n, t) 
and b{n, t). We can reasonable expect that 

a{n, t) + b{n, t) 
2 

will remain close to f{n/N,t), where f{x,t) is a solution to 

In fact, this method should be accurate to order e^. Further, we can expect a whole 
series of conserved quantities similar to (22). 

There are serious problems in discretizing the equations (23) and (24) in time 
efficiently. This problem can be seen by setting a{n) = e^v{n/N) and b{n) = 
e^w{n/N). Then to order e, we have 
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and 

(28) ^=2^-2^. 

dx dx dx 

So while arbitrary solutions tend to move with speeds on the order of magnitude 
1, the special solutions we have constructed with (25) and (26) tend to move with 
speed e^. Thus if one applies some explicit method such as Runge-Kutte, one is 
virtually guaranteed that the differencing in time will produce instability if the 
time step size is not small with respect to 1. One can intuitively understand the 
situation as follows: the equations (23) and (24) define a vector field on the space V 
of functions a and b. For most of the a and b coming from smooth data, the vectors 
in this field have order of magnitude one. But for the set W of those coming from 
(25) and (26), the order of the magnitude of these vectors is e^. Thus W is almost 
fixed by the fiow. In fact, W tends to act as a neutral fixed set, so solutions nears 
W tend to stay near W. A large step size will tend to push away from W, and as 
more steps are taken, the distance from W increases exponentially. In practice, the 
high frequencies are most increased and so instability results. I believe that this 
problem can be overcome using the Crank-Nicholson method. Given a and b we 
wish to find a" and b" , the results of propagating a and b using (23) and (24) over 
a time step At. We try to solve the equations 

a{k) + ^{N^La,b{k)+NFa,b{k)) = a'{k) - ^{N^La^ ,v{k) + NFa' ,v{k)) 

and 

b{k) + ^{N^Ma,b{k) + NGa,b{k)) = b\k) - ^{N^Ma',b'{k) + NGa',b'{k)). 

to obtain an Crank- Nicholson approximations a' and b' to a" and b". If we take the 
step size small with respect of A/"^, then 

is small. So in solving for a' and b' in terms of a and b, the non-linear terms are 
small and so one should be able to find an algorithm to solve these equations using 
some version of Newton's method. At any rate, I tried a rather ad hoc method 
which was stable on the examples I tried. 
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Appendix 

The calculations in this paper were made using the computer algebra system 
Maple V, Release 2. Here is an explanation of the functions in this program, 
followed by the program: 

• dt is a procedure designed to differentiate g by t, given that 

dt 

So dt acts as a derivation of sums, products and exponentials. Further, 

d d 
dt dx 

is replaced by 

d d 
dx dt 

• d(i,n) is di{n). It implements the recursive relation defining the di{n). 

• translate produces the Taylor expansion of f{x + ne) to order N. 

• R is defined in the main part of this paper. 

• aa implements the recursive definition of the Toda equations. 

• XZ[p] is L>i,fe and YZ[p] is L>2,fe. 

• AA[p] and XX are the derivative of / for that particular Toda equation. 

• check[p] checks that the time derivative of B is the same as the derivative 
calculated using the p^^ Toda equation. 

• The last three lines calculate the various linear combinations of AA[p] to get 
the KdV hierarchy. 



dt := proc(g,f ,h) 
local i,n,p; 

if type(g,'+') then 
P := 0; 

n := nops(g) ; 

for i to n do p := p+dt (op(i ,g) ,f ,h) od; 
RETURN (p) 
elif typeCg,'*') then 
P := 0; 

for i to nops(g) do p := p+dt(op(i,g) ,f ,h)/op(i,g) od; 

RETURN (expand (g*p) ) 
elif typeCg,'**') then 

RETURN(op(2,g)*op(l,g)**(op(2,g)-l)*dt(op(l,g),f ,h)) 
elif typeCg, function) and op(0,g) = f then RETURN (h) 
elif type (g, function) and op(0,g) = diff then 
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elif type (g, numeric) or g = epi then RETURN (0) 
fi 

end ; 

d : = 

proc(i ,n) 
options remember; 

if i < then RETURN (0) 

elif i = then RETURN (1) 
elif n = then RETURN (0) 

elif n < then RETURN(d(i,n+l)-A(n)*d(i-l,n)-B(n)*d(i-2,n-l)) 

elif < n then RETURN(d(i,n-l)+A(n-l)*d(i-l,n-l)+B(n-l)*d(i-2,n-2)) 

fi 

end 

trainslate : = 
proc(f ,n) 
local i,y; 

y := f(x); for i to N do y := y+dif f (f (x) ,x $ i)*(epi*n)**i/i ! od 
end ; 

R : = 
proc(x) 

-l/4*diff (f (x) ,x)-l/8*epi*f (x)**2+l/192*diff (f (x) ,x $ 3)*epi**2+ 

epi**3*(l/64*f (x)*diff (f (x) ,x $ 2)+l/64*dif f (f (x) ,x)**2-l/32*f (x)**3)+ 
epi**4*(-l/7680*diff (f (x) ,x $ 5)+l/64*f (x) **2*dif f (f (x) ,x) )+epi**5*( 
3/256*diff (f (x) ,x) **2*f (x)+3/512*f (x) **2*dif f (f (x) ,x $ 2)-5/512* 
f (x)**4-l/i536*f (x)*diff (f (x) ,x $ 4)-3/2048*diff (f (x) ,x $ 2)**2-l/384 
*diff (f (x),x)*diff (f (x),x $ 3)) 

end 

N := 11; 

F := proc(x) 2+epi**2*f (x)-epi**3*R(x) end; 
G := proc(x) -l+epi**2*f (x)+epi**3*R(x) end; 

XZ[1] := B(0)-B(1); 
YZ[1] := B(0)*(A(0)-A(-1)); 
for k from 2 to 4 do 
aa := proc(p,n) 
local r,X; 

X := d(k-p,n+k)-d(k-p,n) ; 

for r from p+1 to k-1 do X := X-aa(r,n)*d(r-p,n+r) od; 
X 

end ; 
XZ[k] := expand(aa(l,0)*B(l)-aa(l,-l)*B(0)) ; 
YZ[k] := expand(B(0)*(aa(0,0)-aa(0,-l))) 

od; 

A := proc(n) options remember; translate (F,n) end; 
B := proc(n) options remember; translate (G,n) end; 
for p to 4 do 

XX := convert (taylorC 1/2* (XZ [p] +YZ [p] )/epi**3, epi, 12) .polynom) ; 
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YY := taylor(dt(G(x) ,f ,XX)-YZ[p]/epi,epi,9) ; 
YY := collect (expand (YY) ,epi) ; 
check [p] := YY; 
print (p , check [p] ) 

od; 

printlevel := 1; 
AA[1] ; 

AA[2] := collect(expand(AA[2]+2*AA[l] ) ,epi) ; 

AA[3] := collect(expand(AA[3]-6*AA[l]+2*AA[2]) ,epi) ; 

AA[4] := collect(expand(AA[4]+20*AA[l]-6*AA[2]+2*AA[3]),epi) 
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